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^ : Abstract 

> 

Q>^ I We derive a new factorization formula in perturbative quantum chromody- 

^D ■ namics for the Drell-Yan massive lepton-pair cross section as a function of 

1^^ ■ the transverse momentum Qt of the pair. When Qt is much larger than 

^D . the pair's invariant mass Q, this factorization formula systematically resums 

the logarithmic contributions of the type a^ In"^ {Qj^ / Q'^) to all orders in the 

^ ^ strong coupling a^. When Qt ~ Q, our formula yields the same Drell-Yan 

Q_i' cross section as conventional fixed order QCD perturbation theory. We show 

Q^' that resummation is important when the collision energy \/S is large enough 

D . and Qt 3> Q, and we argue that perturbative expansions are more stable and 

reliable in terms of the modified factorization formula. 
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I. INTRODUCTION 

The study of massive lepton-pair production in hadronic collisions (the Drell Yan process) 
has been a valuable pursuit for many years 0. The process is an excellent laboratory for 
theoretical and experimental investigations of strong interaction dynamics, and it is a channel 
for discovery of quarkonium states and intermediate vector bosons. In the Drell Yan process, 
the massive lepton-pair is produced via the decay of an intermediate virtual photon, 7*. 
Within the context of perturbative quantum chromodynamics (QCD), the Drell- Yan cross 
section in a collision between hadrons A and B, A{Pa) + B{Pb) -^ 7*(~^ ^KQ)) + ^? can 
be expressed in terms of the cross section for production of an unpolarized virtual photon 
of the same invariant mass , 



dQ^ dQl dy ~ \3nQy dQ^dy 



(1) 



The variables Q, Qt-, and y are the invariant mass, transverse momentum, and rapidity of the 
pair. Symbol X stands for an inclusive sum over final states that recoil against the virtual 
photon. An integration has been performed over the angular distribution in the lepton-pair 
rest frame. Because the leptons can be detected and measured without restrictions, massive 
lepton-pair production as well as inclusive virtual photon production defined in Eq. (|l|) are 
entirely inclusive. 

Precise knowledge of the gluon parton distribution in nucleons is critical for reliable 
predictions of the signals and backgrounds for many important reactions studied at the Fer- 
milab Tevatron and CERN Large Hadron Collider (LHC). It was pointed out recently that 
the transverse momentum distribution of massive lepton-pairs produced in hadronic colli- 
sions is an advantageous source of constraints on the gluon distribution 0, free from the 
experimental and theoretical complications of photon isolation that beset studies of prompt 
photon production [0,^. Other than the difference between a virtual and a real photon, 
the Drell- Yan process and prompt photon production share the same partonic subprocesses. 
Similar to prompt photon production, the lowest-order virtual photon "Compton" subpro- 
cess: g + q ^>- '~i* + q dominates the Qt distribution when Qt > Q/2, and the next-to-leading 
order contributions preserve the fact that the Qt distributions are dominated by gluon ini- 
tiated partonic subprocesses 0. 

If both physically measured quantities Q and Qt are large, the cross section for lepton 
pairs of invariant mass Q and transverse momentum Qt can be factored systematically in 
QCD perturbation theory and expressed as [B,B| 



Y^ I dXi4>a/A{Xl, fJ^) I dX2(f)b/B{x2, /i) 1^^^^^^^'^^' ^' '^^' ^' ^^ ' ^^^^ 



"^^^QldT ^Vj ^^^'^'^Z^^^^' ^) / ^^^^^^^(^2, /i) ^^%?f^^ ' 



The sum J2a,b runs over all parton flavors; (pa/A and (pb/B are normal parton distributions; 
and /i is the renormalization and the factorization scale. The function d(yab-^'y*(Q)x/dQ%dy in 
Eq. d^) represents the short- distance physics of the collision and is calculable perturbatively 
in terms of a power series in Q;s(/i). The leading order and next-to-leading order contributions 
are available @J^. The scale /i is of the order of the energy exchange in the reaction. 



/^ ~ VQ' + Ql- 



There is a phase space penalty associated with the finite mass of the virtual photon, and 
the Drell-Yan factor aem/(3vrQ^) < 10~^/Q^ renders the production rates for massive lepton 
pairs small at large values of Q and Qt- In order to enhance the Drell-Yan cross section 
while keeping the dominance of the gluon initiated subprocesses, it is useful to study lepton 
pairs with low invariant mass and relatively large transverse momenta 0. With the large 
transverse momentum Qt setting the hard scale of the collision, the invariant mass of the 
virtual photon Q can be small, as long as the process can be identified experimentally, and 
the numerical value Q 2> Aqcd- For example, the cross section for Drell-Yan production was 
measured by the CERN UAl Collaboration |^ for virtual photon mass Q G [2r?7,^, 2.5] GeV. 

When Qt is very different from Q while both are much larger than Aqcd? the calculation 
of massive lepton-pair production becomes a two-scale problem in QCD perturbation theory. 
The corresponding short-distance partonic parts, calculated in conventional fixed-order QCD 
perturbation theory, include potentially large terms proportional to the logarithm of the ratio 
of these two physical scales. As a result, the higher-order corrections in powers of a^ are 
not necessarily small. The ratio a^^^ /a^^ [oc a^ x (large logarithms)] can be of order 1, 
and the convergence of the conventional perturbative expansion in powers of a^ is possibly 
impaired. 

When Q^ <^ Q^, the Drell-Yan (or W^ and Z) transverse momentum distributions 
calculated in fixed-order QCD perturbation theory are known not to be reliable 0. After all- 
orders resummation of the large \n{Q'^ /Q^) terms is performed, predictions for the transverse 
momentum distributions become consistent with data for Q\, <^ Q^ |[TO| , p!T| . Similarly, 
when Q'^ ^ Q^, the region of phase space of interest in this paper, the perturbatively 
calculated short-distance partonic parts, dd'ab^'y* (q)x / dQTdy in Eq. (^, receive one power 
of the logarithm \n{Q'^/Q'^) at every order of a^ beyond the leading order. At sufficiently 
large Qt, the coefficients of the perturbative expansion in a^ will have large logarithmic 
terms, and these high order corrections may not be small. In order to derive reliable QCD 
predictions, resummation of the logarithmic terms ln'"(Q|./Q^) must be considered. 

The purpose of this paper is to modify the factorization formula in Eq. (|]) so that 
resummation of the logarithmic contributions can be included naturally when Q^ ^ Q^. At 
the same time, the modified factorization formula should remain effectively the same as the 
conventional factorization formula in Eq. @ when Q^ ~ Q^. 

In the next section, we review the general structure of the Drell-Yan cross section, and 
we identify where the large logarithms arise when Q^ ^ Q^ ^ Aq(-.j-). We show that when 
Q^ is fixed and Q^/Qt ~^ 0, the Drell-Yan cross section behaves similarly to the cross 
section for prompt real photon production 0. The large logarithmic contributions to the 
Drell-Yan cross section at high Qt come from partonic processes that fall into the two-stage 
generic pattern of fragmentation contributions: (1) short-distance production of a parton of 
momentum Pc, and (2) fragmentation of this parton into the observed virtual photon. 

In Sec. |Tl|, we show that the ln"^{QT/Q'^) logarithmic contributions to the Drell-Yan 
cross section can be resummed systematically to all orders in a^. We demonstrate that 
these resummed logarithmic contributions have the same factored form as those for single 
hadron production at large transverse momentum, with the parton-to-hadron fragmentation 
functions replaced by the fragmentation functions for a parton to a virtual photon of invari- 
ant mass Q. As for single hadron production, the short-distance production of the parton is 
evaluated at a single hard scale (~ Qt), and it is calculable perturbatively in a power series 



of Us- However, unlike the parton-to-hadron fragmentation functions, the parton-to- virtual- 
photon fragmentation functions are calculable perturbatively if Q^ ^ ^qcd \]2l 



In addition to the logarithmic contributions, the Drell-Yan cross section includes large 
non-logarithmic contributions, in particular, the leading order contributions, referred to of- 
ten as direct contributions. In Sec. |IV|, we derive our modified factorization formula for 
massive lepton-pair production, Eq. (|35D , in which both logarithmic and non-logarithmic 
contributions are fully included. This modified factorization formula agrees with the con- 
ventional fixed-order QCD expression when Q\ ~ Q^ (or when the logarithmic contributions 
are less important). We show that at the next-to-leading order, the difference between the 
modified factorization formula and the conventional factorization formula is completely de- 
termined by QCD evolution of the virtual photon fragmentation functions. Our modified 
factorization formalism reorganizes the single perturbative expansion of conventional QCD 
factorization into two perturbative expansions plus the perturbatively calculated parton-to- 
virtual photon fragmentation functions. The main advantage of this reorganization is that 
the new perturbative expansions are evaluated at a single hard scale and are free of large 
logarithmic terms for Qt > Q- 

In Sec. 0, we present our predictions for the cross sections for massive lepton-pair produc- 
tion at energies of interest for experiments at the Fermilab Tevatron, Brookhaven's RHIC, 
and the CERN Large Hadron Collider. We include both leading order and next-to-leading 
order direct short-distance contributions and the resummed logarithmic contributions. The 
resummed large logarithmic contributions change the shape of the predicted Qr-spectrum of 
the Drell-Yan cross section but, at the order in perturbation theory at which we work, they 
have only a modest effect on the normalization. We confirm that after the large logarithmic 
terms are resummed to all orders in a^ the Drell-Yan cross section at large Qt remains an 
excellent source of contraints on the gluon parton density. Our conclusions are summarized 



in Sec. VL 



II. MASSIVE LEPTON PAIR PRODUCTION AT FIXED ORDER 

In hadronic collisions, massive lepton pair production proceeds through partonic hard- 
scattering processes involving initial-state quarks and gluons. If the lepton pair's invariant 
mass Q and its transverse momentum Qt are both much larger than Aqcd, the partonic hard- 
scattering at a distance scale between 0{1/Q) and 0{1/Qt) can be systematically factored 
from the physics at the scale of hadron wave functions, 0(1/Aqcd)- In this situation, the 
cross section can be expressed in the factored form of Eq. (^. Corrections to the expression in 
Eq. (|^) are suppressed by powers of Aqqj-,/(5^ or K^^q-q/Qt- The predictive power of Eq. (^ 
relies on the universality of the parton distributions and the reliability of the partonic cross 
sections. 

The short- distance partonic cross sections, daab^Y iQ)^ / ^Qt^V in Eq- (0), can be calcu- 
lated in principle order by order in QCD perturbation theory in a power series of the strong 
coupling as, 

dQ\dy =]:^Hl,L^^xi^„x,,Q,QT,y;f.) [-^j . (3) 



The reliability of QCD perturbative calculations depends on the behavior of the coefficient 
functions H^_^^,x i^i Eq. (|^). 

At lowest-order, 0{a^g), the only partonic subprocess for virtual photon production is 
g + g — i> 7*. The incoming partons are assumed to be collinear to their respective incoming 
hadrons if power suppressed corrections are neglected. Therefore, the lowest order coefficient 
function, H^^_^^*x ^ ^{Qt), vanishes if Qt 7^ 0. 

At order 0{as)-, both quark-antiquark annihilation, g + g — > 7* + (7, and "Compton", 
(7 + g ^ 7* + g, subprocesses contribute to the Drell-Yan cross section, with the recoil of the 
final-state parton balancing the transverse momentum of the lepton pair. These partonic 
subprocesses, known as the 2^2 subprocesses, are shown in Figs. |I|(a) and (b). They are 
often referred as the leading order (LO) contributions to the Drell-Yan cross section with 
finite transverse momentum. The corresponding leading order coefficient functions are 



^(1) ^ 2 ^ 



9s 



(Svr^) 6 {{x^Pa + x,Ps - Qf) (^) . (4) 



The incoming parton flavors "a6" can be either gg for the quark-antiquark annihilation or 
gq for the Compton subprocess, and e^ is quark's fractional charge. In Eq. @, g^ and agm 
are the strong coupling constant and the flne structure constant of QED, respectively; Pa 
and Pb are the momenta of the colliding hadrons, and 5* = {Pa + Pb)^ is the square of the 
total collision energy. The expressions {Mab^^'xl"^ in Eq- ® are the squares of the matrix 
elements of the partonic subprocesses shown in Fig. 0, summed (averaged) over the colors 
and spins of the flnal-state (initial-state) partons. They are available in Ref. [|T^]. As long 
as Qt is large, the LO coefficient functions in Eq. (^) are well-behaved, even when Q^ — » 0. 
The calculation of the perturbative coefficient functions at order O(a^), known as the 
next-to-leading order (NLO) contribution, involve all 2 — >■ 3 partonic subprocesses with the 
virtual photon in the final-state as well as the 2 — > 2 diagrams in Fig. |l] with one-loop 
corrections. After renormalization, the loop momentum integrations for the 2^2 dia- 
grams at order O(a^) yield renormalization scale (/x) dependence and logarithmic terms in 
the coefficient functions H^f^ *x- Integration over the phase space of the extra parton in 
the final-state of the 2 — i> 3 subprocesses leads to a collinear divergence when this parton 
is collinear to either incoming parton. QCD factorization and subtraction of the collinear 

divergence results in factorization scale (/i/) dependence and logarithmic terms in the co- 

(2) 
efficient functions at this order. Consequently, the coefficient functions H^^i^'^^x display 

logarithmic dependence on the ratios of the following momentum scales: /i, /xj, Qt-, and Q 

m- 

Since we are interested in identifying the logarithms of the ratio Q'r/Q^-, with Qj- ^ Q^, 
we concentrate on the part of the Drell-Yan cross section that diverges when \n{Q^/Q'^) — >■ 
00 with Q^ ^ ^QCD- According to the QCD factorization theorem, the perturbatively 
calculated partonic cross sections daab^'y* {q)x / (^Qt^V in Eq. (||) should be analytic functions 
of Qt and Q^. Therefore, we expect the logarithmic behavior of the Drell-Yan cross section 
as \ia{Q^/Q^) -^ cxd with Q^ ^ ^qcd fixed to be connected closely to logarithmic divergences 
associated with the massless photon (Q^ = 0) in the case of prompt real photon production. 

Other than the non- vanishing invariant mass, production of a virtual photon and a real 
photon share the same parton-level Feynman diagrams. However, the QCD factorization 
formula for production of a real photon (Q^ = 0) is different from that in Eq. (^, 



daAB- 



nx 



dQ^dy 



^ / dXi(f)a/AiXi, H) / (iX206/B(x2,yU) 
a,b 



dQ\ dy dQ'^ dy 



(5) 



where da^^jZ^x / '^Qt'^V represents the direct production of the real photon at a short-distance 
scale of 0{1/Qt)] and 



da. 



dQ^dy 



,,-^ r dz da, 



*cX 



dp^ 



ct 



dy 



Pc 



Q 2 



Dr 



*yX 



>,/i|) 



(6) 



is the fragmentation contribution to the prompt photon cross section. The function 
d&ab^cx / dpl^dy in Eq. (H) represents short-distance production of a parton of flavor c, 
and Dc^^x {z, fi'p) is a fragmentation function for parton c to fragment into a real photon 
with photon momentum Q = zp^, Hf is the fragmentation scale. 

The fragmentation contribution arises because there are collinear singularities associated 
with the region of phase space in which the real photon is parallel to one or more of the 
final-state partons [|13|,|14[. Because of the real photon is massless, the parent parton, which 
fragments into the real photon and other collinear partons, can propagate for a long time. 
Consequently, quantum interference between the production of the parent parton and the 
physics associated with the fragmentation (or decay) of the parton is suppressed. Therefore, 
the fragmentation contribution to prompt photon production can be further factored as in 
Eq. (P). Because the transverse momentum Q^ is large, all logarithmic collinear divergences 
associated with the massless photon arise from final-state partons that are parallel to the 
observed real photon. Such logarithmic divergences are all absorbed into the fragmentation 
functions. These functions are nonperturbative in nature. 

Unlike prompt photon production, it is not necessary to introduce fragmentation func- 
tions to absorb final-state collinear singularities for the Drell-Yan cross section. Because the 
photon is off-shell, its large invariant mass Q regulates the singularity. This finite mass reg- 
ularization leads to a logarithmic dependence of the Drell-Yan cross section on the invariant 
mass of the virtual photon Q. If Q is large enough and Qt, the only other physically ob- 
served momentum scale, is not too large, the logarithmic terms ln™'(Q^/Q^) are small, and 
no resummation of the logarithms is necessary for a reliable prediction of the cross section. 

When Q^ is chosen to be small, so as to enhance the Drell-Yan cross section, and when the 
collision energy yS and Qt become large, it is necessary to examine the size of the final-state 
logarithmic contributions and ascertain whether resummation of these logarithmic terms is 
warranted. 

The explicit form of NLO contributions to the Drell-Yan cross section from a 2 ^ 3 
partonic subprocess g + g — >• 7*(Q) + q' + q' [§] provides an example of the logarithmic 
terms under discussion. The relevant Feynman diagrams are shown in Fig. ^(a) and 0(b) for 
initial- and final-state photon radiation. Since we are interested mainly in the structure of 
the final state in the region where a photon becomes collinear to a quark, this subprocess is 
typical of the generic 2^3 subprocess in Fig. ^, and conclusions drawn from it can applied 
to other subprocesses such as qg and qq' scattering. For the phase space integrals of the 
diagrams in Fig. Q, we use the exact forms in Ref. [0. The logarithmic contributions have 
the form 

s + Q^ -S2 + X' 



In 



s + g2 _ S2 - A 



g{s,t,u,Q' 



The function g{s, t, u, Q'^) is given in Ref. 0. It is well-behaved as Q^ — > 0. The parton-level 
Mandelstam variables s, t, u are defined as s = {pi +^2)^, i = {'Pi~ QYi and u = {p2 — QY- 

Q? 



The function A = J {t + uY — 4(^^52, with S2 = (pi + P2 — QY = s + t + u — Q^\ S2 is the 
square of the invariant mass of the two final-state partons that recoil against Q. 

In the limit Q^ <^ |t + m|, the generic logarithm from the splitting of a quark into a 
photon can be approximated as 



In 



S + Q^ -S2 + \ 



In 



2Q' 



Q')/{-{t + u)) 



In 



2(-(t + «)) 

-{t + uY 



Q' 



(7) 



To derive the second line we use Q^ -C |t + u\ and s ~ |t -|- m|. Since |t + «! ~ 0{Q^ + Q"^), 
the logarithm takes the form ln[(Q| + Q'^)/Q'^] ~ ln[Q|/(5^], and the limit of Q"^ <^ \t + u\ 
is effectively the same as Q^ -C Q^- The logarithmic contributions to the Drell-Yan cross 
section become more important when Q^ is chosen to be small. 

In principle, we can resum these large but finite logarithms into parton-to- virtual-photon 
fragmentation functions, just as the logarithmic divergences in real photon production are re- 
summed into real-photon fragmentation functions. Because the logarithms in the Drell-Yan 
case are finite, and Q^ is much larger than Aq(jj-), the parton-to- virtual-photon fragmentation 
functions should be calculable perturbatively |T^ . 



III. FRAGMENTATION CONTRIBUTIONS TO THE DRELL-YAN CROSS 

SECTION 

In this section, we derive the resummed logarithmic contributions to the Drell-Yan cross 
section when Q^ ^ Q"^. Since we are interested mainly in the cross section at low Q^, we 
ignore contributions from the intermediate vector boson Z. 

A. Logarithmic Contributions to the Drell-Yan Cross Section at Large Qt 

As demonstrated in the last section, the Drell-Yan cross section at large transverse 
momentum Qt receives potentially large logarithmic terms ln'"(Q|./Q^) from the part of 
phase space in which the virtual photon is almost coUinear to one or more final-state partons. 
Since at least one final-state parton is needed to balance the virtual photon's transverse 
momentum, the logarithmic contributions can arise only at NLO and beyond. 

Because of the logarithms, the coefficient functions in Eq. (j^) might be large, and resum- 
mation of the logarithmic contributions might be needed. As long as Q^ is much larger than 
Aq(jj-), all coefficient functions in Eq. (^) are calculable in principle order by order in QCD 
perturbation theory. Therefore, resummation of the logarithms In^lQ^/Q'^) is actually a 
reorganization of the perturbative expansion in Eq. (^, such that all coefficient functions in 
the reorganized perturbative expansions are evaluated at a single hard scale and free of any 
large logarithms. 



The energy exchange in the hard colhsion is of the order of JQj^ + Q^ ~ Qt + 0{Q'^ /Q\). 
When Qj' ^ Q^i the partonic hard colhsion should not be sensitive to the scale Q^ at which 
the virtual photon is produced. If we neglect power corrections of the order Q^ IQ\, the 
partonic hard collisions are effectively independent of Q^, except for the logarithmic depen- 
dence ln™'((5y/Q^) from the final-state bremsstrahlung production of the virtual photon. 
Therefore, other than the appearance of the virtual photon's mass Q to regulate the final- 
state coUinear divergences, the potentially large logarithmic contributions at large Qt have 
the same structure as the fragmentation contributions to prompt real photon production. 
They can be separated into two stages: (1) production of a parton of momentum pc at a 
very short- distance (~ 1/]3ct ~ 1/Qt), and (2) production of the lepton-pair via a virtual 
photon of invariant mass Q through bremsstrahlung (or fragmentation) from the parton pro- 
duced at the first stage. This two-stage production, shown in Fig. § or in general in Fig. ^, 
shares the same generic pattern of the fragmentation production of a single particle (e.g., 
a hadron of mass M^, or a real photon) at large transverse momentum Q^. If we neglect 
the power suppressed quantum interference between these two stages, the fragmentation (or 
bremsstrahlung) contributions should have the same general factored form that is present 
in single hadron or prompt photon production 



dQ\dy 



E 



dz 



da. 



ab—*cX 

dpl^ dy 



Xi,X2, Pc 



Q/z; jj,f) 



Dr 



.r*x 



^z, fip] Q ). 



(8) 



Superscript (F) indicates the fragmentation contribution; the sum J2c runs over all parton 
flavors; and ftp is the fragmentation scale defined below. The four vector Q^ in Eq. (|]) is 
defined to be Q'^ but with Q^ set to be zero. 

The partonic cross sections, ddab^cx / dpl^ dy in Eq. (^, represent the inclusive produc- 
tion of a parton of flavor c. These partonic cross sections are evaluated at a single hard scale 
PcT ~ Qt- Exact perturbative expressions for the partonic cross sections depend on how 
the fragmentation functions Dc-,.y*x{z, fip] Q^) are defined (or depend on the choices of the 
scheme) [|12| . 

As for single hadron (or prompt photon) production, all large logarithms from the inte- 
gration over the distance scale from 0{l/pcj,) (or 0{1/Qt)) to 0{1/Q), where the virtual 
photon is produced, are resummed into the fragmentation functions. The main difference 
for production of a virtual photon, a hadron, or a real photon is the difference in the frag- 
mentation functions. 



B. Parton-to- Virtual-Photon Fragmentation Functions 



The fragmentation function for a parton of flavor c into a virtual photon of invariant 
mass Q is defined as a probability density for finding the virtual photon with fraction z 
of the parent parton's momentum. Because of the universality of fragmentation functions, 
virtual photon fragmentation functions can be derived in a process independent way |T2[ or 
extracted from a specific physical process. In this subsection, we derive the virtual photon 
fragmentation functions by resumming the leading logarithmic contributions to the Drell- 
Yan cross section at large transverse momentum. 



In order to resum the large logarithmic contributions to the Drell-Yan cross section to all 
orders and to extract the virtual photon fragmentation functions, we must identify all sources 
and the pattern of these logarithmic contributions. Consider a non-singlet fragmentation 
of a quark into a virtual photon, shown in Fig. ^ The virtual photon is produced from 
bremsstrahlung of a quark (or an antiquark). The quark (or antiquark) itself is produced 
either in the hard collision, as shown in Fig. ^(a), or from the fragmentation of another 
quark, as shown in Fig. |^(b) or Fig. |§(c). In a physical gauge (such as the light-cone gauge), 
the leading large logarithmic contributions of the subprocesses in Fig. |^ come from the part 
of phase space in which the daughter quark's invariant mass kf is much smaller than that 
of the parent quark k'^_^_l. In this situation, we can neglect all contributions suppressed by 
powers of kf/kf_^_i. The "decay" rate (or the fragmentation function) of the parent quark 
into a daughter quark plus massless partons is dominated by logarithmic contributions that 
are proportional to ln{kf_^_i/kf). The leading logarithmic contributions from the general 
fragmentation diagram in Fig. ^ arise from the region of strong ordering in the invariant 
masses of the fragmentation partons in Fig. |^: p^ ^ ■■■kf+i '^ k"^ '^ ■■■ '^ k^ ||T2 . 



With the strong ordering approximation, the leading logarithmic contributions from the 
non-singlet fragmentation in Fig. ^ can be factored into the hard production of a quark of 
momentum p^ convoluted with a sum of all-orders fragmentation ladder diagrams, as shown 
in Fig. ^ By comparing the factored cross section in Eq. (^ with the factored expression in 
Fig. p, and summing the ladder diagrams to all orders, one obtains the leading contributions 



to the non-singlet quark-to-virtual-photon fragmentation functions [|12 






>-?^«')^/lt(l7^-^--^»^«^') 



yz k^ 

" r'^hi dkf fa, r^ dzi 




zi,kl;Q')y (9) 



The superscript 'WS"' represents the non-singlet contribution; the upper limit of integra- 
tion is fc^_,_^ = p^; and the lower limit of integration, Q'^/z, is the mass threshold (or mini- 
mum invariant mass) for the quark to produce the virtual photon of invariant mass Q ^2 . 
The leading-order quark-to-quark splitting function P^^ Jz) is the same as the leading- 
order quark-to-quark splitting function of the DGLAP evolution equations |l^ . In deriving 



Eq. (|^), we include the diagrams with quark wave function renormalization in addition to 
the ladder diagrams shown in Fig. ^j. The function 'yq_^^*{z, fcg; Q^) in Eq. (H) is the lowest 
order QED splitting function for a quark to fragment into a virtual photon 0,|l3]; defined 
below. 

Since the invariant mass of the parent quark, p^, can be very large in high energy colli- 
sions, the resummed logarithmic contributions, given by the second term on the right-hand- 
side of Eq. d^), can be important for the Drell-Yan cross section at low Q^. From Eq. (H) and 
the fact that the partonic hard parts are evaluated at a single hard scale ~ pcj, ~ Qt, we 
conclude that all leading large logarithmic contributions at high Qx are included in the vir- 
tual photon fragmentation functions. Resummation of the large logarithmic contributions is 
equivalent to the derivation of the virtual photon fragmentation functions in Eq. (pi). Unlike 



the parton-to-hadron fragmentation functions, the virtual photon fragmentation functions 
in Eq. (|^) have no dependence on any non-perturbative momentum scale. The virtual pho- 
ton's non- vanishing invariant mass removes the final-state collinear singularities that appear 
in the parton-to-hadron (or real photon) fragmentation functions. Therefore, all parton-to- 
virtual-photon fragmentation functions should be free of collinear singularities, and they are 



calculable in principle perturbatively to all orders in a^ [|12 |. 

By re-organizing the second term on the right- hand- side of Eq. (^), one may derive an 
integral equation for the non-singlet quark-to-virtual-photon fragmentation function. 



^^ (y 'r?-n'^\ - r^ ^^0 /"em (0) f^ ;^2. ^2^ 



fPc dk"^ ( a., [^ dz' „,n^ . z 

2/2 



-/;;^ferT-.'24).-..(.',*^;Q,^ (10, 



Carrying out the sum over all-orders logarithmic contributions in Eq. (^ is equivalent to 
solving the integral equation in Eq. (|T0|). 

We introduce fragmentation scale iip and let /xl^ = vli the square of the invariant mass 
of the parent quark in Eq. (p!oD . One can derive an evolution equation by applying fipd/dfip 
to both sides of Eq. ([10|), 

+ 5^ r f^n'S,(^)flfirx(-', fh Q') ■ (11) 

Because the quark can interact directly with the virtual photon, the evolution equation in 
Eq. ( pTl) has an inhomogeneous term. 

By extending the simple ladder diagrams in Fig. ^ to general ladder diagrams ITS 



one 



can derive evolution equations for the singlet quark-to-virtual-photon and gluon-to-virtual- 



photon fragmentation functions [|T2 

A 
d\i 



f^F-fT^c~^rx{z,fil]Q'^) = (^y^J7c_T,.(z,/i|,a^;Q^ 



+ (?-) E /' ^^c-..(4, «s) D,^rx{z\ /il; Q') , (12) 



.27r. .... 

where c,d = q,q,g- In Eq. (|I^) , the evolution kernels Pc^d are evaluated at a single hard 
scale, fiF, and can be calculated perturbatively as a power series in a^. QCD corrections 
to the QED quark-to-virtual-photon splitting function '~fc-,'y* can be evaluated in principle 
order-by-order in a^. 

Calculating the lowest order quark-to- virtual-photon ladder diagram in Fig. ^, we obtain 
the leading order quark-to- virtual-photon QED evolution kernel [|^,0, 



i%{z,e-Q')=el 



1 + (1-^)^ JQ' 
z \ zk'^ 



e{e - ^) . (13) 



The ^-function is a consequence of the mass threshold. The gluon-to-virtual-photon evolu- 
tion kernel vanishes at the same order 
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j^^X,,{z,e;Q'] 







(14) 



because the gluon does not interact directly with the virtual photon. 

It is important to note that if we work in the strict leading power (or leading twist) 
approximation, we would drop both the power corrections to the fragmentation functions as 
well as power corrections to the evolution kernels of the fragmentation functions. That is, 
we would neglect the 0{Q^/k'^) term in Eq. (p!3D , 



7, 



(LP-O) 



{z,e) 






k'; Q' = 0)= el 



e{e- 



Q' 



(15) 



The superscript "LP" represents the leading power approximation. With this strictly leading 
power quark-to-photon QED evolution kernel, the evolution equation in Eq. (|l^) for the 
parton-to- virtual-photon fragmentation functions is exactly the same as that for the parton- 
to- real-photon fragmentation functions [l^. Under the strict leading power approximation, 
the only difference between virtual and real photon fragmentation functions is the boundary 
condition for the evolution equations. For real photon fragmentation functions, a set of 
unknown non-perturbative input fragmentation functions Dc-^^x{z) is needed at a given 
scale yUp, (~ a few GeV) [|19[. On the other hand, no non-perturbative input fragmentation 
functions are needed to solve the evolution equations in Eq. (|T^. Instead, the mass threshold 
for production of a time-like virtual photon of invariant mass Q imposes a natural boundary 
condition for all flavors c, 



Dr 



^7*X 



>,/x|<QV^;g2) = 0. 



(16) 



The strict leading-power approximation might be too severe in the threshold region [12 



Because its mass is non-zero, the virtual photon can have both transverse and longitudinal 
polarization modes. The QED evolution kernel in Eq. (^3|) is a sum of evolution kernels for a 



quark to fragment into either transverse (T) or longitudinal (L) polarization modes ||T2|,|T7 



with 



itAz, k'; Q') = 2 7f4(., e-, Q') + i^:i,. {z, e-, Q\ 



JO) /. u2.r^2^ 



(0) 



,Jz,k';Q' 
z,e;Q') 



1 /1 + (1 



zk^ 



9{k' 



z 



Q 



Q' 



Jkr^^^^'-V- 



(17) 

(18) 
(19) 



The factor 2 in Eq. (p!7| ) represents the two transverse polarization states of the virtual 
photon. Under the strict leading power approximation, Eq. ([I9|) vanishes. In this strict 
limit, only transversely polarized virtual photons are produced through the fragmentation 
processes, and there are no logarithmic contributions to the production of longitudinally 
polarized virtual photons. Furthermore, without the 0{Q'^/k'^) term in Eq. ([T8|) , the evo- 
lution kernel 7q~i^* i^^ ^^i Q'^) gives the wrong threshold behavior. Instead of being zero at 
the threshold when k"^ = Q'^/z, the leading power evolution kernel 7^^^* {z, k"^; Q"^) is finite 
and proportional to 1/ z. It is large if z is small. 
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In general, there can be two types of power corrections to the leading power virtual 
photon fragmentation functions. Power suppressed contributions 0{Q'^/fi'j^) to the frag- 
mentation functions are one type (Type-I as defined in Ref. [0), and power corrections 
contributions to the evolution kernels (or the slopes) of the fragmentation functions are the 
other type (known as Type-II). For example, consider the lowest order contribution to the 
quark-to-virtual-photon fragmentation function from the evolution kernel in Eq. (Ill 










In 



Q"^ ) \ zjjj, 



(20) 



The (1 — Q"^ / z^y) term results from the power suppressed Q^ j zk^ term in the evolution 
kernel, and it is clear that this term is as important as the logarithmic term in the threshold 
region. The term proportional to 1 in this {1 — Q'^/zfi'^) combination is not power suppressed 
by 0{Q'^ / ii\). In the threshold region, the virtual photon fragmentation functions are 
dominated by the longitudinally polarized component particularly when z is small |]12||. The 



Type-II power corrections are not necessarily small and could provide power non- suppressed 
contributions to physical observables fill . 



In the rest of our discussion, we keep the leading power suppressed terms in the QED 
evolution kernels in Eq. (0) when we calculate our parton-to-virtual-photon fragmentation 
functions. As shown in Ref. |T^, the inhomogeneous QED evolution kernels in Eq. (|T2p 



dominate the scale dependence of the fragmentation functions, and therefore, we neglect 
the power corrections to the QCD evolution kernels Pc^d iii Eq. (|T2p. With the inclusion 
of power corrections in the evolution kernels, the resummation discussed here is no longer 
a simple one-scale problem in QCD perturbation theory. More detailed discussions of the 
virtual-photon fragmentation functions can be found in Ref. |1^ . 



In summary, the all-orders resummation of the large logarithmic contributions to the 
low mass Drell-Yan cross section is equivalent to a sum of all logarithmic contributions to 
the virtual photon fragmentation functions, achieved by solution of the evolution equations 
in Eq. (|T^. The evolution equations for the virtual photon fragmentation functions have 



the same functional forms as those for real photon fragmentation functions |jT3|. However, 
the differences in both boundary conditions and the inhomogeneous terms due to the non- 
vanishing of Q^ lead to many differences between the real and virtual photon fragmentation 
functions. One major difference is that the virtual photon fragmentation functions are purely 
perturbative. In addition, the virtual photon fragmentation functions provide the resummed 
contributions to the production of a longitudinally polarized virtual photon |12| . 



C. Calculation of the Partonic Hard Parts 

To determine the fragmentation contributions to the Drell-Yan cross section. 



/7 (^) 



a,b 



dQ^dy 
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X 



E 



dz 



dpl^ dy 



Xi,X2, Pc = Q/z]fiF 



Dc^-y*x{z, fJ'Fl Q 



(21) 



we must evaluate the partonic hard parts da^^l^^-^ / dpl^dy in Eq. (p]) . Although these 
are calculable perturbatively, exact expressions depend on how the parton distributions 
and fragmentation functions are defined (or the choice of factorization scheme). In this 
subsection, we provide a self-consistent procedure for the calculation of the partonic hard 
parts. 

We separate the procedure into four steps: (1) instead of considering the hadronic pro- 
cess: A + B ^ 7*(Q)X, we apply the factored formula in Eq. (^) to a partonic pro- 
cess: a' + b' —^ c'X with an on-shell final-state parton c'; (2) we expand both sides of the 
factored formula for the partonic process order-by-order in a^; (3) we calculate the par- 
tonic cross section daa'b'-^c'x/dpli dy, and parton-to-parton distributions and fragmentation 
functions order-by-order in a^; and (4) we extract the short-distance partonic hard parts 
d&ab^cx/dpl^dy by comparing both sides of the perturbatively expanded factored formula 
at the same order of a^. 

By applying the fragmentation expression in Eq. ( |2lD to the partonic processes: a' + b' — »• 
c'X, we obtain the following schematic formula 



0"a'b' 



Z^ <Pa/a' 



%/h' 



^ab-*c 



Dr 



(22) 



Symbol C?> represents the convolutions over momentum fractions Xi, X2, and z in Eq. (pTI). 
(In Eq. (P^p and in the rest of the equations of this subsection, we omit the inclusive symbol 
X for reasons of notational simplicity.) 

To produce a parton with large transverse momentum, a 2 — > 2 partonic process is 
required at the parton level, at minimum of the order O(a^). Expanding both sides of 
the factored formula Eq. (^2]) order-by-order in as, we define the following perturbative 
expansions. 



C^a'b' 



a 



(F) 
ab^c 



V- (n) (ajjiy\ 

n=2 V 27r ; 

/ y a 



ab—->c 



n=2 



n=0 \ 



\ 27r 
asif^y 



D, 



n=0 \ 



277 ; ' 

(n) /as(/^)\" 



2tt 



We substitute these four perturbative expansions into Eq. (^2]) and obtain at 0{a 






a 



.(2) _ V- , (0) 

a'b'^c' — Z^ y^a/a' 



^b/b' 



H. 



ab^c 



D 



(0) 



(23) 
(24) 
(25) 
(26) 

(27) 



Since the zeroth order parton distributions and fragmentation functions are 5-functions, 



Eq. 



yields Hf^j_^^ 



a 



(2) 
ab— 



or 
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dpl^ dy dpl^ dy 



(28) 



at leading order. 

Expanding Eq. p^ ) to NLO, we write 



a 



'b'^c' — 7 , Va/a' 


® 4/i. ® Hfl 


.®fl™. 


+ EC 


® 4/1. ® ^i'f- 


o^ijfi 


+ y:e 


» ■*!,;!,. ® ffS^ 


.■SDfl 


a,b,c 


» <!>',% ® ffS^ 


.^Dfl 



(29) 

Using the zero-th order parton distributions and fragmentation functions, for 

dd-[^Zc^'^^/dpl^ dy, we obtain 

o-(3) _ ^(3) V- , (1) rr{2) 

a' 

- E <p'^jb ® H^Lc - E ^Slc ® ^?lc- (30) 



r(2) 



In Eq. (|30|) , H^^Uc is the leading order contribution calculated in Eq. (^81) . The partonic cross 

(3) 

section (Tab_^i. and all parton-level parton distributions and fragmentation functions are per- 



turbatively calculable with proper regulators. The subtraction terms in Eq. ( pOD remove the 
collinear singularities associated with the massless partons. Following the same procedure, 
we can derive the short- distance partonic hard parts for the fragmentation contributions in 
Eq. ( p]) at all orders in a^. 

Equation ( |30D shows that beyond the leading order, the exact expressions for the pertur- 
batively calculated short- distance hard parts depend on the definitions of parton-level parton 
distributions and fragmentation functions. The partonic hard parts are fixed uniquely once 
we fix the parton-level parton distributions and fragmentation functions. In order to use 
available conventional parton distributions, we have little choice other than to select the 



parton-level parton distributions in either the MS or DIS scheme ||20 

In the partonic cross section cr^^Lc^ the parton momentum pc is assumed to be massless, 
pI = 0, and, therefore, a final-state collinear singularity arises. Within the usual QCD 
factorization framework, there is freedom to choose any factorization scheme to remove 
the final-state collinear singularities of the partonic cross section and absorb all possible 
finite differences into the non-perturbative fragmentation functions. Different choices for 
the factorization scheme lead to finite differences between the extracted non-perturbative 
fragmentation functions. However, owing to the non-zero invariant mass of the virtual pho- 
ton, the parton-to-virtual-photon fragmentation functions do not have final-state collinear 
singularities. They are completely perturbative. Therefore, the parton-to-virtual-photon 
fragmentation functions cannot uniquely fix the definition of the parton-to-parton fragmen- 
tation functions. As a consequence of the difference in the invariant masses of the parton 
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and the virtual photon, an extra constraint has to be introduced to specify the parton-to- 
parton fragmentation functions D]!!'\^c- ^^ ^^^ choose a scheme for defining the parton-level 
fragmentation functions so as to remove the final-state collinear singularities in the par- 
tonic cross sections cr^l^ci with n > 3. The finite differences between schemes cannot be 
completely absorbed into the perturbative parton-to- virtual-photon fragmentation functions. 
Finite differences associated with the choice of scheme for the parton-level fragmentation 
functions correspond to non-logarithmic contributions to the Drell-Yan cross section. For the 
logarithmic contributions to the Drell-Yan fragmentation functions, we can choose parton- 
level fragmentation functions in the MS scheme or in any other scheme when calculating 
the partonic short-distance hard parts in Eq. (pQ]). The perturbatively calculated partonic 
hard parts da^^i^^^/dpl^dy will be the same as those for inclusive single hadron (or prompt 
photon) production if the same factorization scheme is used. The non-logarithmic differences 
caused by the different choices of the factorization schemes can be absorbed into the direct 
contributions to the Drell-Yan cross section, defined and discussed in next section. 



IV. DRELL-YAN CROSS SECTION WITH RESUMMED FRAGMENTATION 

CONTRIBUTIONS 

In this section, we derive our modified factorization formula for the Drell-Yan cross 
section including resummation of the logarithmic contributions. 

Particularly at leading order, but also at all higher orders in perturbation theory, there 
are significant contributions to the cross section that are not included in the fragmentation 
terms. It is essential that the full factorization formula include all non-logarithmic contri- 
butions order-by-order in a^, while retaining the all orders resummation of the logarithmic 
contributions. 

We use the word direct to designate the non-logarithmic contributions to the Drell-Yan 
cross section . This use of the term is the same as its use in prompt photon production. At 
leading order the direct terms in the Drell Yan cross section are those supplied by the short- 
distance 0{1/Qt) 2-parton to 2-parton Compton and annihilation subprocesses. However, 
there is an additional component at higher orders. The direct contribution must also absorb 
the finite non-logarithmic differences in the partonic hard parts Hc^l,cx associated with the 
ambiguity in definition of the parton-level fragmentation functions, as discussed in the last 
section. This second component accounts for the non-logarithmic terms over the distance 
interval 0{1/Q) to 0{1/Qt)- The logarithmic terms are included in the fragmentation 
functions. The physics of our direct term is very similar to that represented by the Y term 
in the CSS formalism [Q for resummation at small Qt- 

To be precise, we define the direct contribution as the difference 



dQj^ dy dQj- dy dQj- dy 

= ^ / dXi(j)a/Aixi, H) / dX2(pb/Bix2,fl) 



a,b 






dQj^dy 
The parton level direct term is obtained from the factored formulas in Eqs. (0) and (^ 
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(31) 



, - (Dir) 



{xi,x2,Q,QT,y;fJ',fJ'F) 



da, 



with da, 
tional 



(F) 



/dQ'^dy given in Eq. 
QCD 



fixed-order 



afe->7*(Q)X 

dQj'dy 

(F) 
afc^7*(Q)X 

dQ\dy 
and da, 



da 



{xi,x2,Q,QT,y;fJ') 



xi,x2,Q,QT,y;fJ',fJ'F) 



(32) 



.(F) 



ab-*-t*{Q)xldQ\dy calculated in conven- 
perturbation theory. Since 



daah-.'y*{Q)xldQ'^dy and da^f^l^^^Q.-^/dQ^dy are calculable perturbatively, the direct con- 



tributions, da[^Zr^Q)x 



/dQ'ip dy in Eq. 
have the perturbative expansion 



should also be calculable perturbatively. They 



, - {Dir) 
"'^ah^^'{Q)X 

dQ'^dy 



n=\ 



(n) 
afe^7*(Q)X 



[xl,x2,Q,QT,y■,^J'F,^J') 



2tt 



(33) 



where n > 1 because the leading order contributions to Drell-Yan cross section are 
of order 0{aemOis)- Since the fragmentation contributions da\^J ,,Q^-^/dQ\dy and 
daah-^-y* {Q)x / dQ"^ dy share the same large logarithmic terms, ln™(Q^/Q^), the direct contri- 
butions should be free of large logarithms order- by-order in a^. Therefore, the coefficient 
functions ^^^ .(qw in Eq. ( |55| ) are evaluated at a single hard scale (~ Qt)-, and the per- 
turbative expansion for the direct contribution in Eq. (|33D is expected to be well-behaved 
perturbatively. We substitute Eqs. (§, (El), (E§, and (0) into Eq. (P^) and obtain 



Y, 



(n) 



ah~^-^*(Q)X 



H. 



(n) 



ab^j*(Q)X 



E 

m=2 



E^< 



ab-^cX '^ ^ 



r{Q)x 



(34) 



{n—m) 



Functions D)!'''y,Q^■^ are the coefficient functions for the perturbatively calculated parton-to- 
virtual-photon fragmentation functions |l^. The subtraction term in Eq. (|34D removes the 
logarithmic contributions included in the coefficient functions -ffab^7*(Q)x that are calculated 
in conventional fixed-order perturbation theory. As remarked above, the direct contributions 
defined in Eq. (B3) absorb the finite non-logarithmic differences in the partonic hard parts 



H. 



(m) 
ab- 



cX- If-f^ab 



(m) 



^cX 



is larger in one scheme and more is included in the fragmentation contri- 



butions, then more will be subtracted from the direct contributions. Any non-compensated 
differences will be higher order corrections in as- 

We substitute Eq. ( ^2]) into Eq. (|3T|) and derive a modified factorization formula for the 
Drell-Yan cross section at Qt > Q'- 

daAB^'y*{Q)X 



dQ^dy 



E 

a,b 



dXi(l)a/AiXi, fl) / dX2(l)b/Bix2,fi) 



,^{Dir) 
"'^ab^r{Q)X 




{xi,x2,Q,QT,y;fJ'F,fJ') 



(35) 



(F) 
ab^rcX , 



[Xi,X2, Pc 



Both perturbatively 

[Dir) 
'afe^7*(Q)X/ 



calculated 



CT^y 

short- distance 



z 



flF,fJ') 



Dr 



'{Q)x{z,^p;Q 



hard parts. 



d^^ablcx/dp% dy 



and da.i .„ir^\v/dQ'^dy, are free of large logarithms. All potentially large \n!^ {Q\ / Q"^) 
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terms in the region Q^ ^ Q^ are resummed into the virtual photon fragmentation functions 
Dc^^*{Q)x- The direct contribution (first term on the right-hand side of Eq. (|35D) represents 
the production of lepton-pairs at the distance scale of the hard collisions (~ IJQt)- The 
fragmentation contribution stands for the sum of all leading logarithmic contributions from 
a distance scale at \/Qt to l/Q. 

We note here, in passing, that unlike the logarithmic terms in the region of small Qt, 
the logarithmic terms at large Qt are typical DGLAP logarithms associated with collinear 
contributions. At low Qt-, logarithmic contributions arise from the collinear and the infrared 
regions. Correspondingly, in resummation formalisms relevant at small Qt-, there are leading 
log terms (including both collinear and infrared logs at all orders in a^), and next-to-leading 
log terms (including only one of the two types of logs), and "next-to-next-to- leading" log 
terms, and so forth. In our case, we do not have the type of "next-to-leading" logs that 
appear in resummation procedures at small Qt- Anything left-over in the large Qt region, 
after resummation of the leading logs of DGLAP type, is included in the direct terms. 

Since we have in mind applications at large Qt but small Q, it is important to consider 
the possible role of higher- twist contributions proportional to inverse-powers of Q. Like 
power corrections from target mass effects, power corrections in our case should appear 
in the form of m^/ (energy exchange)^, with m ~ Q. Since Q'^ sets the hard scale, and 
Q^ acts as a mass threshold, the only dimensionless ratios we expect to see are Q^ /Q\ 
and Q^ I ^ 1 but no l/Q^ term. There can be three kinds of power corrections: power 
corrections to the partonic hard parts, suppressed by Q^ /Qt] power corrections to the 
fragmentation contributions, suppressed by Q'^/fij? ~ Q'^/Qt'i ^^^ power corrections to the 
evolution kernels of the fragmentation functions, proportional to Q'^/fi^. Since /i^ runs from 
Q^ to fijp, the third type could be significant, Q^//i^j„ ~ 0(1). Although this third variety 
is potentially important, the terms should be much smaller than the leading logarithmic 
contributions as long as ln(Q^/Q^) ^ 1. 

Our modified factorization formula for the Drell-Yan cross section in Eq. ( P5| ) is very 
similar to that for prompt real photon production in Eq. (^. However, the Q^-dependence 
in the direct production term and differences in the fragmentation functions distinguish the 
Drell-Yan virtual photon production from prompt real photon production. 

The key difference between our modified factorization formula in Eq. (|35D and the con- 
ventional factorization formula in Eq. (0) resides in the way the logarithmic contributions 
from final-state parton splitting are handled. Instead of one perturbative series in powers 
of as in the conventional approach, we have two perturbative expansions in our modified 
factorization formula: one for the direct and one for the fragmentation contribution. All co- 
efficient functions in the new perturbative expansions are free of large logarithms. The large 
logarithms in the conventional perturbative expansion are resummed into the fragmentation 
functions. In the rest of this section, we demonstrate the difference by comparing the two 
factorization formulas explicitly order- by-order in powers of ««. 

According to the conventional factorization approach, the lowest order partonic con- 
tributions to virtual photon production are provided by the quark-gluon Compton and 
quark- antiquark annihilation diagrams at order aemCts^ 
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7 (C-LO) 
^'^AB^-y*{Q)X 



J2 J dXi(l)a/AiXi, H) J C?X206/B(x2,/i) 



,^(C-LO) 

dQj:dy 



{xi,x2,Q,QT,y;fJ') 



(36) 



The superscript [C — LO) stands for the leading order conventional perturbative calculation. 
In our modified factorization approach, Eq. (^), there are two perturbative expansions 
corresponding to the direct and fragmentation contributions, respectively. Each perturbation 
series has its own leading order contributions. Because the subtraction term in Eq. (^) starts 
at 0(af), the LO direct contribution is the same as the LO conventional contribution in 
Eq. (HD, 



, (Dir-LO) 

dQ'^dy 



E 

a,b 



dXi(j)a/A{xi, fi) / dx2(j)b/B{x2, fJ-) 



,^ {Dir-LO) 

dQ'^dy 



{xi,x2,Q,QT,y;fJ') 



(37) 



The LO partonic hard part da\f^^^,Q<^x/dQTdy equals t/a 



(C-LO) 

ab^-y* {Q)X/ 

The virtual photon fragmentation functions are of order Oem 



/dQ^dy in Eq. (|6]). 
and the partonic hard 

parts, da^fjl^^x / dPcT ^H ^^ ^^- dH) start at order a^. The LO fragmentation contributions 
to the Drell-Yan cross section would seem therefore to be of order 0(aem«s); one power of 
as higher than the LO direct contributions in Eq. (|H7|). However, the large logarithms from 
the virtual photon fragmentation functions are proportional to ln(/i|,) oc l/asifip), meaning 
that the LO fragmentation contributions to the Drell-Yan cross section can be viewed as 
effectively 0{aem(^s) terms, of the same order as the LO direct contribution: 



, (F-LO) 
^^AB^"I*{Q)X 

dQ^ dy 



Y^ / dXi(j)a/AiXi, IJ,) / dX2(t)b/B{.X2,lj) 



^^ fdz da[lJ^^x^ 



dpL dy 



xi,X2, Pc= —; Hf, ^) Dc^y,(Q)x{z, uliQ^) . (38) 



In our modified factorization formalism, the LO contribution to the Drell-Yan cross section 
at large Qt is equal to the sum of the LO terms of the two perturbative expansions: 



da 



(LO) 

AB^T'{Q)X 

dQj^dy 



jADir-LO) (jJP-LO) 

"''-' AB^-i*{Q)X ""AB-»7*(Q)X 



dQj'dy 



dQ^dy 



(39) 



The LO direct and fragmentation contributions are found in Eq. 

tively. 



and Eq. 



respec- 



At NLO, aem«s5 i^ ^^e conventional fixed-order approach, the Drell-Yan cross section has 



the same form as that in Eq. (|36D but with superscripts {C — LO) replaced by {C — NLO). 
According to our modified factorization formalism, the NLO term should be the sum of 
the NLO terms in both the direct and fragmentation contributions in Eq. (^5l). The NLO 
direct contribution is 



, (Dir-NLO) 
"^AB-*7*(Q)X 

dQj^ dy 



J2 J dXi(j)a/AiXi, IJ,) J (iX20b/ij(x2,/i) 



a,b 



X 



,^{C-NLO) 

ab^y {Q)x ^3,^^ ^2, Q, Qt, y; /i) 



dQ^ dy 

E 



(40) 



dz da^^-^°^ 
z^ dpl^ dy 



{Xi,X2, Pc 



Q 



;/^f,/^)-D^^^*(q)jy(^5/^f5Q ) 



The subtraction term is a consequence of the definition of the direct contribution in Eq. (p4D. 
The subtraction term is also necessary to remove the logarithmic terms from the conventional 
NLO expression in Eq. (^OD, to avoid double counting. These logarithmic contributions 
have been included in the LO fragmentation contribution to the Drell-Yan cross section in 
Eq. (H). 

We would need NLO fragmentation contributions to complete the calculation of the NLO 
contribution, 



, (F-NLO) 

dQ\dy 



Y^ I dXi(f)a/A{xi, fi) I dX24>b/B{x2, /i) 

r dz da^^-^^°^ 6 



a,b 



dPcr- dy 



/iF, Ai) Dc^^*(^Q)xiz, fip] Q^ 



(41) 



The virtual-photon fragmentation functions include the all orders resummation of the log- 
arithmic terms, and they are the same as those for the LO contribution in Eq. (0). The 
0(a|) partonic hard parts dal^_^^x /dplj, dy are defined in Sec. |IIIC| and will be presented 
elsewhere. 

The NLO contribution to the Drell-Yan cross section at large transverse momentum is 



, (NLO) 
^^AB->-i*(Q)X 

dQ'^dy 



, (Dir-NLO) , {F~NLO) 



dQ^dy 



dQf-dy 



(42) 
and CETI), 



where the NLO direct and fragmentation contributions are found in Eqs. (| 
respectively. 

In our modified factorization formalism, the LO and NLO contributions are different 
from those in the conventional formalism because the virtual photon fragmentation functions 
include all orders of the large leading logarithmic contributions. The difference is better seen 
if we rewrite our LO and NLO results in the following form. 



da 



AB^7*(Q)X 



/7 (-f^O) 
"''^AB^j*{Q)X 



, (NLO) 



dQj^ dy dQj, dy dQj^ dy 

= E / dXi(j)a/A{Xl, fJ^) I dX2(f)b/B{x2, /i) 



(43) 



a,b 



X 



da 



(C-LO) 
ab~^j*{Q)X 

dQj. dy 



{xi,X2,Q,QT,y]fJ') + 



,^{C-NLO) 
^ ^ab^r{Q)X 

dQj'dy 



xi,x2,Q,QT,y;fJ') 
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a,b 

r r dz da^^-^'^^ O 



z^ dpl^ dy - - z 



X "^ , , 

, (F-NLO) 



r dz da^^^^""^^ Q 1 



On the right-hand-side of Eq. (^41), the first term the prediction of the conventional fixed- 
order factorization formahsni in Eq. (Q); the second term is the difference between our 
modified factorization formula and the conventional fixed-order factorization formula up to 
NLO. 

To conclude this section, we emphasize that our modified factorization formalism 
in Eq. (^) effectively reorganizes the single perturbative expansion of conventional 
QCD factorization, daah^'y*(Q)x/dQ'^dy, in Eq. (^, into two perturbative expansions, 

dd'abZi*(Q)x/'^QTdy and dali,_^^x/dplj,dy, plus the perturbatively calculated evolution ker- 
nels. The main advantage of this reorganization is that the new perturbative expansions 
are evaluated at a single hard scale and are free of large logarithms. As shown in the next 
section, the ratios of the NLO over the LO contributions in the new perturbative expansions 
are smaller than the ratios evaluated in the conventional approach. 



V. NUMERICAL RESULTS AND PREDICTIONS 

In this section, we present numerical evaluations of the leading and next-to-leading order 
Drell-Yan cross section at large transverse momentum. We show the quantitative differences 
between the predictions of the modified factorization formula of Eq. (^) and the conven- 
tional factorization formula. We demonstrate the sensitivity of the cross section to the gluon 
distribution at low Q^ and high Qt- 

We employ the CTEQ5M set of parton distributions [^. We use a two-loop expression 
for the strong coupling strength a^, with the value of Aqcd specified by CTEQ5M, and a 
one loop expression for aem with aem{Mz) = 1/128. We choose Q^ as the renormalization 
scale for aem- We equate the renormalization and factorization scales and set hard scale 



fi = fif = kJQt + Q"^ with constant n = 0(1). 

The normal MS factorization scheme removes the ultraviolet 1/e pole of the parton- 
level fragmentation functions along with the corresponding splitting functions for massless 
partons, but the scheme does not guarantee that the fragmentation functions to a massive 



parton or to a photon with no n- vanishing invariant mass [|T^] will be positive. As a result, the 
virtual photon fragmentation functions calculated in the MS scheme can be negative in the 
region of large z [p!2|,p!7|. As long as the cross section is positive, a negative fragmentation 



function is simply a particular separation of finite contributions between the coefficient 
function and the fragmentation function. It is not a problem in principle. Nevertheless, it 
is more appealing intuitively that the fragmentation functions be positive definite. Since 
the virtual photon fragmentation functions are purely perturbative, it is possible to preserve 
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positivity of the fragmentation functions if we require that the mass threshold constraints 
be respected at every stage of the fragmentation (or bremsstrahlung radiation). 

As shown in Ref. |]12| , the ultraviolet (UV) divergences of the virtual photon fragmen- 



tation functions come from two sources: (1) elementary divergent diagrams associated with 
the renormalization of the fields and coupling constants in QCD and QED, and (2) the loop 
momenta of the skeleton ladder diagrams. An invariant mass cutoff scheme is introduced 
in Ref. |T2[ to render the fragmentation functions positive definite. In this scheme, all UV 



divergences associated with the internal elementary divergent diagrams are removed in the 
same way as the QCD and QED Lagrangian are renormalized (say, in the MS scheme). 
The UV divergences connected with loop momenta of the skeleton ladder diagrams and 
wave function renormalization of the composite operators are removed by imposition of an 
invariant mass cut on the loop momenta. All running coupling constants in this scheme are 
renormalized in the same way as the renormalization of the Lagrangian. Use of the same 
invariant mass cutoff on the loop momenta of the skeleton ladder diagrams and the virtual 
diagrams due to the wave function renormalization of the composite operators ensures the 



infrared cancellation between the real and the virtual diagrams |T6|. In this invariant mass 
cut-off scheme, a parton-to- virtual-photon fragmentation function Dc-^^*(q)x{z, fip; Q^) can 
be viewed as an inclusive rate for "decay" of the parton of flavor c and squared invariant 



mass fXp into a virtual photon of squared invariant mass Q and momentum fraction z [jT2 

In our numerical calculations, we use the virtual photon fragmentation functions in the 
invariant mass cut-off scheme of Ref. [|12[ . The invariant mass for quark c in Fig. |^ to decay 



into a virtual photon of invariant mass Q and a massless quark is 

— * 

The three-vector Qt is perpendicular to the direction of the parent quark's momentum. 
The fragmentation scale is chosen to be the invariant mass of the fragmenting parton. For 
example, we choose /i|, = p^ at lowest order. Equation ( ^5] ) shows that the mass threshold 
requires that the fragmentation scale be ^p > Q'^/z |lT2| , p!7| . If we keep the perturbative con- 
tributions to a high enough order, the cross section should not be sensitive to the choice of 
fragmentation scale. However, since the fragmentation scale in this new scheme is very dif- 
ferent from the traditional scale in the MS scheme, we test two choices for the fragmentation 
scale. 



fXF = K ^JQI + Q2 ^ (46) 

/zp = «:^M±^. (47) 

The first choice is the same as that for the renormalization and factorization scales. 

The second choice in Eq. (^7p is motivated by the fact that the squared invariant mass 
of the fragmenting parton is of the order {Q\. + Q"^)/ z, obtained as follows. For the generic 
2 — > 3 diagram in Fig. ^, with the final parton that recoils against Pc assumed massless, we 
compute 

pI = 2{pa + Pb) ■ Pc - (Pa + Pbf 

?a . (48) 
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To derive the second line, we use the approximations Q^ <^ Q^ and s < \t + u\/z. If we let 
pI, Eq. (^) leads to the fragmentation logarithm, 



/^l~ 



In 






In 






In 



\t 



u 



Q' 



(49) 



u\ ~ 



0{Ql + Q'), we find fil 



consistent with the logarithm in Eq. (|^. Since \t 
{Ql + Q^)/z. 

The perturbatively calculated partonic parts in conventional perturbation theory, 
^ab^-i*(Q)x ^^"^ ^ah^'y*(Q)x ^^ -^*1- PI)? ^-^e available and calculated in the MS scheme P,0. 
The fragmentation contribution in Eq. ( pS]) depends on the leading order short- distance 
partonic parts, 0"^5^cX ? ^^^ on the virtual photon fragmentation functions. The partonic 
parts are 



.{F-LO) 



TT 



dpi dy 2X1X2S 



-^Mab^cX 



8n^) 5 ({xiPa + X2PB - Pcf 



2tt 



(50) 



The functions Mab^cx are the lowest order matrix elements for partons a and h to produce 
a parton of flavor c, averaged over the colors and spins of the partons in the initial state. 
They are available in Ref. 



g. 



To compare with data, we introduce the invariant differential cross section 



E 



d^CJAB- 



>i+i-X 



d^Q 



TC 



a^ 



dQ^ 
dQ 



da 



AB^e+e-(Q)x 
dQ^ dQl dy 



37r2 



Q' 



da 



AB^-f*(Q}X 

dQj'dy 



(51) 



The integration in dQ"^ is over a bin centered on the invariant mass Q. The differential cross 
sections daAB^i+i-x/dQ^ dQ^dy and daAB^-y*x/dQ'^dy are given in Eq. (P and Eq. (^3|), 
respectively. If the bin size AQ is much smaller than Q, the cross section can be approxi- 
mated as 



E 



d^aAB- 



>e+e-x 



d^Q 



2aem AQ 
37r2 Q 



da 



AB^j*{Q)X 

dQ\dy 



(52) 



The differential cross section through next-to-leading order for production of a virtual pho- 
ton, daAB^'y*{Q)x/dQTdy, is given in Eq. (|3]). 

Since we do not present the NLO partonic hard parts for the fragmentation contribu- 
tions in this paper, we use Eq. (^3]) without the da^^^^^ /dpl^dy term for the numerical 
calculations of the Drell-Yan cross sections. 

In Fig. 0, ^ and ^, we plot the differential cross section Ed^a /d^Q as a function of Qt 
for pp — > £+^^ + X at the Tevatron collider energy v^ = 2.0 TeV, for pp — ^ £"*"£" + X at 
the LHC energy yS = 14 TeV, and for pp -^ i^£~ + X at the RHIC proton-proton energy 
y/S = 500 GeV. In all three figures, the rapidity y = 0. The renormalization, factorization, 
and fragmentation scales are set equal, with constant k = 1. The solid, dashed, dotted, 
and dot-dashed lines represent the total, the LO direct, the NLO direct, and the resummed 
fragmentation contributions, respectively. We show results for virtual photon invariant mass 
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2 < Q < 3 GeV and 4: < Q < 5 GeV. (Note: the solid curve represents the sum of the 
dashed, dotted, and dot-dashed hues.) To enhance the cross sections, the smaller value of 
the invariant mass Q is favored as long as the lepton pairs can be identified experimentally. 
At the integrated luminosities of run I of the Fermilab collider, prompt real photons have 
been observed with values of transverse momentum extending to 100 GeV and beyond [^ 




Scahng to the massive lepton pair case, we judge that it should be possible to examine cross 
sections in the same data sample out to Qt of 30 GeV or more for virtual photon invariant 
mass 2 < Q < 3 GeV. Values of Qt to 50 GeV or so may be reached with 2fb~^ at run II. 
To demonstrate the renormalization, factorization and fragmentation scale-dependence 
of the cross sections, we introduce the ratio 



(53) 



The denominator in Eq. ([53| ) is obtained with all three scales equal: fi = fij = ftp = 

JQt + Q^- III Fig- HO' ^^ plot -R/x as a function of Qt at the Tevatron energy and y = 0, 
for 2 GeV and 5 GeV, with different scale choices. We first fix all three scales to be the 
same: fi = fif = fip = ^^JQt + Q^i with k = 2 (solid) and k = 1/2 (dashed). The dotted 

lines correspond to the choice /i = /i/ = JQt + Q^, and ^p = v (Qr + Q^)/^; as defined in 
Eq. (^). As shown in Fig. |T0|, the scale dependence is not great. It yields an uncertainty of 
about 15 % for all reasonable values of Qt at the Tevatron energy. The effect of the different 
choice for the fragmentation scale in Eq. (^7]) is also small. The same features are preserved 
at the other collider energies. 

To show the quantitative difference between the conventional fixed-order perturbative 
calculations and the calculations with all-orders resummation of the large logarithmic terms, 
we introduce the ratios 



RiLO) ^ 



, (Dir-LO) 

d^Q 



d^Q 



d^Q 



(54) 



for the LO contributions, and 



R 



(Dir-LO) 



(Dir-NLO) 



(F-LO) 



jp '^" AB^l+l-X _|_ jp "'" AB^i+l-X _|_ j^"'" AB-^l+l-X 



d^Q 



d^Q 



d^Q 



jp "'^ AB^l+l-X _|_ jp^_^ AB^i+t-X 



d^Q 



d^Q 



(55) 



The ratio R^^'^^ is the ratio of the LO contributions in the two different factorization 
formalisms. Since the LO direct contribution is the same as the LO term in the conven- 
tional calculation, any deviation of R^^^^ from unity measures the size of the logarithmic 
contributions. In Fig. 11, we plot the ratio R^^^' as a function of Qt at different energies 
for 2 < Q < 3 GeV and 4 < Q < 5 GeV. The solid, dashed, and dotted lines represent the 
Tevatron, LHC, RHIC energies. Owing to the threshold behavior of the fragmentation func- 
tion, the fragmentation contribution vanishes for Qt ~ Q-, and R^^^^ -^ 1. The logarithmic 
contributions are very important at the LHC energy and less important at RHIC energies. 
The reason for this energy dependence is the phase space penalty associated with the large 
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invariant mass of the virtual photon. The large logarithm in the fragmentation function is 
proportional to \iii{z fijp / Q"^) for each power of a^. Since the strong coupling strength ^^(/i) 
is proportional to l/ln(/i^/AQQj-,), the product 

«.(/.) ln(z/.|/Q^)ocMfii^ (56) 



becomes of order unity only when fip ^ fi^ Q^. Otherwise, the combination in Eq. (|56|) is 
not large because of the factors z and Q^ from the mass threshold. 

The ratio R in Eq. (^) measures the difference between the conventional calculation up 
to NLO accuracy and our resummed calculation, obtained from our modified factorization 
formalism without the NLO fragmentation contributions. Without the NLO fragmentation 
contributions Eda^^~^^'^^ / cP'Q in the numerator, the ratio R does not represent the entire 
ratio of the NLO contributions in the two different factorization approaches. Nevertheless, 
this ratio does indicate some of the effects of QCD resummation. 

The difference between the numerator and the denominator of the ratio R in Eq. ( |55D 
is determined entirely by the second term in Eq. (|44D without the NLO fragmentation 
contributions, 

f f f dz da^^~^'^^ Q 

^J J ' J z^ dpi^dy z 

X [D,^,.^Q)x{z,ixl;Q')-D^X'^Q^xi^,ixl;Q')] . (57) 

This difference is proportional to the difference between the LO QED parton-to-virtual- 

'"^ ( 



photon fragmentation functions D\,]^^,q.^{z, ^\]Q'^) and the corresponding QCD evolved 



parton-to-virtual-photon fragmentation functions Dc-^y*(^Q)x{z, fip;Q ). As shown in 



Ref. |lj], one of the major differences between these fragmentation functions is the behavior 



at large z: 

D,^^*(Q)x(z^l,/i|;Q') = 0. (58) 

QCD evolution reduces the fragmentation function at large z while it increases the fragmen- 
tation function at small z. If the cross section is dominated by the small (large) z region, 
the numerator of the ratio R in Eq. (^) is larger (smaller) than the denominator calculated 
in fixed-order perturbation theory. 

The z-integration of the Drell-Yan cross section runs from Zmjn to 1. We introduce a 
cutoff Zc to limit the integration to the range z^i^ to Zc, and we define the ratio 

R. ^ e'^^^^^{z^) E^^^^^iz^ = 1) . (59) 

The cutoff Zc can be between z^in and 1; Rz{zc = ^^min) = 0, and Rz{zc = !) = !. The shape 
of the ratio Rz establishes which region of z dominates the z-integration. 



In Fig. 12, we plot Rz as a function of Zc at the Tevatron energy with Q = 2 GeV. 



The solid and dashed lines correspond to transverse momenta Qt = 5 and 50 GeV. If the 
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integrand for the ^-integration were independent of 2;, the ratio Rz would be proportional 
to {zc — 2min)/(l " z^in) i which corresponds to a straight line for R-^. The shape of R-^ in 
Fig. |T2| shows that the 2;-integration is dominated by the large z region. 

Since the z-dependence of the partonic cross section in Eq. (l38| ) tends to cancel the 1/ z^ 
factor in the 2;-integration, the 2;-integration for the Drell-Yan cross section is determined 
mainly by the convolution of parton distributions and fragmentation functions and by the 
shapes of these functions. 

For a given collision energy vS, the transverse momentum Qt = z p^j, (x zx y/S with 
X ~ Xi ~ a;2 for y ~ 0. For a fixed value of Qt, a large (small) z corresponds to the 
small (large) x region. Because the parton distributions and fragmentation functions are 
both steeply falling functions of the momentum fractions {x or z), and the Drell-Yan cross 
section is proportional to two powers of parton distributions and one power of fragmentation 
function, the convolution of parton distributions and fragmentation functions favors the 
combination of a small x and large z. The fragmentation contributions are dominated by the 
small X and large z region. The z-integration over the difference between the QCD evolved 
fragmentation function and the leading order QED fragmentation function in Eq. ( |57D is 
dominated by the large z region. The net contribution from this term is expected to be 
negative, and the ratio R in Eq. (|^) should be less than one for Qt ^ Q- 

We plot the ratio R in Fig. O as a function of Qt at y = for 2 < Q < 3 GeV and 



4 < Q < 5 GeV at different collision energies. The solid, dashed, and dotted lines in Fig. ^ 
correspond to the Tevatron energy v^ = 2.0 TeV, the LHC energy v^ = 14 TeV, and 
the RHIC proton-proton energy y/S = 500 GeV, respectively. The ratio R is less than one 
when Qt > Q, as expected. As Qt grows, z^i^ becomes larger, and the fragmentation 
contributions are dominated by even larger values of z. Consequently, i? is a decreasing 
function of the transverse momentum Qt- The predicted distribution in Qt is steepened 
relative to the conventional expectation, with the steepening being greater for smaller \/S. 
The shape of the ratio R represents the effect of QCD evolution on the virtual photon 
fragmentation functions. As Qt — > Q, -R — >• 1, in keeping with the expectation that our 
modified factorization formula should yield the same cross section as conventional fixed-order 
perturbation theory when Qt ~ Q- The apparent slight difference of R from 1 near Qt = Q 
may be attributed to the choice of fragmentation scale fip. If we were to set /i|. = Q%/z, 
our result would be the same as the conventional one when Qt = Q. 

The main advantages of our modified factorization formula in Eq. (|35|) are that the large 
logarithms are removed from the coefficient functions of conventional factorization and all 
perturbatively calculated coefficient functions are evaluated at a single hard scale. The 
reliability of the perturbative calculations is enhanced. In order to demonstrate this benefit, 
we introduce the ratios 

, (C-NLO) I , (C~LO) 

^' = ^ d^Q I ^ d^Q ' ^^°^ 

, (Dir-NLO) I , (Dir-LO) 

TD — jp "'"AB^e+e-x / T? AB^i+i-x ((y^\ 

^^- = ^ d^Q I ^ d^Q ' ^^^^ 

D — zp AB^f+f-x / Tp "''^AB^e+e-x /Ko^ 

^^ = ^ — d^ — / ^ — d^ — • ^^^^ 

The subscripts c, Dir, and F represent the conventional perturbative calculation, our direct 
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contribution, and our fragmentation contribution, respectively. Because the large logarithms 
are removed in the perturbatively calculated coefficient functions in the modified factoriza- 
tion formula, we expect the ratios Roir and Rp to be smaller than the ratio Re- 



We compare Roir and Re in Figs. |14| and |15| at the Tevatron and LHC energies. As 
expected, Roir is smaller than Re. We remind the reader that in our notation, NLO cor- 
responds to the pure O(a^) contribution and LO corresponds to pure 0{as)- Thus, the 
commonly defined i^-factor is obtained by adding 1 to the ratios Figs. |T1| and [T^. The fac- 
tor K would be near 2.3 for Qt ~ 10 GeV (near 1.9 for Qt ~ 50 GeV) at VS = 2 TeV and 
2 < Q <3 GeV in the conventional approach, but somewhat smaller, ~ 2 and ~ 1.6, for our 
direct contribution. The relatively smaller role of the NLO term lends greater confidence in 
the reliability of the calculated Drell-Yan cross sections in our modified formalism. 

The denominators are identical in the definitions of the ratios Re and Roir in Eqs. (|60D 



and (plf ). The difference between the solid and dashed curves in Figs. ^ and |T^ is therefore 
the ratio of the logarithmic piece in NLO to the leading-order direct contribution. The 
figures show that a considerable fraction of the NLO correction in conventional perturbative 
calculations comes from logarithmic contributions associated with bremsstrahlung radia- 
tion of the virtual photon. As noted above, the NLO fragmentation contributions are not 
available yet, and, therefore, we do not show results here for Rp- 



Taken at face value, the ratio R in Fig. 13 suggests that the cross section predicted in our 



modified formalism is modestly less than that in the conventional approach for Qt > Q and 
that the shape of the predicted Qt distribution is steeper. However, we remind the reader 
that we have not included the NLO fragmentation contributions. When these are added 
in subsequent work, the predicted Drell-Yan cross sections from our modified factorization 
formula should be larger than those from the conventional NLO perturbative approach. Es- 
timating the ratio of the fragmentation contributions by the ratio of the direct contributions, 
Rp ~ Roir, we expect the Drell-Yan cross sections, with the NLO fragmentation contribu- 
tions included, will be greater by 10 to 30 % than the predicted totals in Figs. ^ |^- Foi" 
example, at -v/S* = 2 TeV and 2 < Q < 3 GeV, we estimate increases of 15% and 10% at 
Qt = 9 and 39 GeV. 

To demonstrate the relative size of gluon initiated contributions, we define the ratio 

P _ daAB^r {Q)x (gluon-initiated) / daAB^r{Q)x /^o\ 

'~ dQldy I dQldy ' ^ ' 

The numerator includes the contributions from all partonic subprocesses with at least one 
initial-state gluon, and the denominator includes all subprocesses. In Fig. |T^, we show Rg as 
a function of Qt at ?/ = for the Tevatron energy (solid), the LHC energy (dashed), and the 
RHIC proton-proton energy (dotted) at Q = 2 GeV and Q = 5 GeV. We confirm that gluon 
initiated subprocesses dominate the Drell-Yan cross section and that low mass Drell-Yan 
lepton-pair production at large transverse momentum is an excellent source of information 
on the gluon distribution 0. The falloff of Rg at very large Qt is related to the reduction 
of phase space and the fact that cross sections are evaluated at larger values of the partons' 
momentum fractions. In the Tevatron case, quark- ant iquark subprocesses are expected to 
(re) assert their dominance at very large Qt because of the valence nature of the antiquark 
density in the incident antiproton. 
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VI. CONCLUSIONS 

In this paper we introduce a new factorization formalism for the Drell-Yan cross section 
at large Qt- It incorporates all-orders resummation of large logarithmic contributions into 
parton-to-virtual photon fragmentation functions. In addition to the logarithmic contribu- 
tions, our expression for the cross section includes the usual non-logarithmic contributions, 
referred to often as direct contributions, in both leading 0{as) and higher orders. This mod- 
ified factorization formula agrees with the conventional fixed-order QCD expression when 
Qt ~ Q^- The difference between the modified factorization formula and the conventional 
formula is determined by QCD evolution of the virtual photon fragmentation functions. 
Our modification reorganizes the single perturbative expansion of conventional QCD fac- 
torization into two expansions plus the perturbatively calculated parton-to-virtual photon 
fragmentation functions. The new expansions are evaluated at a single hard scale and are 
free of large ln'"(Q^/Q^) terms when Q\. ^ Q^. The ratio of the next-to-leading 0{a1) to 
leading order 0{as) contributions is smaller in the new expansion than in the conventional 
case. The reliability of the predictions is enhanced. 

As shown in Fig. |TT], the contributions from the logarithmic terms are important. They 
can be as large as 70% of the non-logarithmic contributions at the LHC energies. At RHIC 
energies, the logarithmic contributions are smaller because phase space at large transverse 
momentum is limited and the penalty associated with large invariant mass Q is felt more 
strongly. In Fig. |^, |^, and ^, we present predictions through next-to-leading order at en- 
ergies of interest for experiments at the Fermilab Tevatron, Brookhaven's RHIC, and the 
CERN Large Hadron Collider. Resummation somewhat steepens the shape of the predicted 
QT-spectrum relative to the conventional fixed-order approach. However, at the order in 
perturbation theory at which we work, it has only a modest effect on the normalization. 
We confirm that the Drell-Yan cross section at large Qt remains an excellent source of 
constraints on the gluon parton density. 
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FIG. 1. Feynman diagrams for the LO contribution to the Dreh-Yan cross section: (a) 
quark-antiquark annihilation q + q ^ 'J* + g, and (b) Compton g + q ^ ^* + q subprocesses. 
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FIG. 2. Feynman diagrams that ihustrate situations in which the photon can become colhnear 
to a quark in (a) the initial-state and (b) the final-state. 
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FIG. 3. A generic diagram of the lowest order 2^3 subprocess that contributes to the 
Drell-Yan cross section with large final-state logarithmic terms when Q^ <^ Q^. 
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FIG. 4. Sketch of the fragmentation contribution to low mass Drell-Yan lepton-pair production 
at high Qt- 
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FIG. 5. Scattering amplitudes that provide large logarithmic contributions to the Drell-Yan 
cross section via quark fragmentation. 
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FIG. 6. Factored non-singlet quark-to-virtual-photon fragmentation contributions to the 
Drell-Yan cross section. 
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FIG. 7. Drell-Yan cross section as a function of Qt at the Tevatron energy \/5 = 2.0 TeV 
and rapidity y = for the mass intervals (a) 2 < Q < 3 GeV and (b) 4 < Q < 5 GeV. Solid, 
dashed, dotted, and dot-dashed lines stand for the total, LO direct, NLO direct, and resummed 
fragmentation contributions, respectively. 
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FIG. 8. Drell-Yan cross section as a function of Qt at the LHC energy vS^ = 14 TeV and 
rapidity y = for the mass intervals (a) 2 < Q < 3 GeV and (b) 4 < Q < 5 GeV. Sohd, 
dashed, dotted, and dot-dashed hnes stand for the total, LO direct, NLO direct, and resummed 
fragmentation contributions, respectively. 
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FIG. 9. Drell-Yan cross section as a function of Qt at the RHIC proton-proton energy 
^S = 500 GeV and rapidity y = for the mass intervals (a) 2 < Q < 3 GeV and (b) 
4 < Q ^ 5 GeV. Solid, dashed, dotted, and dot-dashed lines stand for the total, LO direct, 
NLO direct, and resummed fragmentation contributions, respectively. 



32 



±.^ 






1.15 


,.-"'" 




1.1 


VS=2 TeV 




1.05 

1 

0.95 


/ Q=2GeV 




- /\ ^°) 




0.9 


'■ \.~___________ 




0.85 


- 




0.8 


,,,,,,,,,,,,,,,,,,,,,, ,,,,,,i,,,,i,,,,,,, 






5 10 15 20 25 30 35 40 45 50 

Qt (GeV) 



±.^ 




1.15 




1.1 


VS=2 TeV 


1.05 
1 


Q=5GeV 


\ 


0.95 


- \ (b) 


0.9 


^^~^^-~--______ 


0.85 


- 


0.8 


,,,,,,,,i,,,,i,,,,i,,,,i,,,,,, ,,,,,,,,, ,,,,,,,, 



5 10 15 20 25 30 35 40 45 50 

Qt (GeV) 



FIG. 10. Ratio R^ in Eq. ( |5^ ) as a function of Qt at y = and the Tevatron energy v5 = 2 TeV 
for virtual photon invariant mass (a) 2 GeV and (b) 5 GeV. The curves are explained in the text. 
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FIG. 11. Ratio w- ' in Eq. (|5^ as a function of Qt at y = and three different collision 
energies for the mass intervals (a) 2 < Q < 3 GeV and (b) 4 < Q < 5 GeV. The curves are 
explained in the text. 
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FIG. 12. Ratio Rz in Eq. ( |59[ ) as a function of Zc at the Tevatron energy \/S = 2 TeV, y = 0, 
and (5 = 2 GeV. The sohd and dashed hnes correspond to transverse momenta Qt = 5 and 50 GeV. 
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FIG. 13. Ratio R in Eq. (^) as a function of Qt at y = and mass intervals (a) 2 < Q < 3 GeV 
and (b) 4 < Q < 5 GeV. Solid, dashed, and dotted lines are for the Tevatron energy y^ = 2.0 TeV, 
the LHC energy v5 = 14 TeV, and the RHIC proton-proton energy v5 = 500 GeV, respectively. 
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FIG. 14. The ratios in Eqs. (60) and ( |6l|) as a function of Qt at y = and for the Tevatron 
energy \/5 = 2.0 TeV and mass intervals (a) 2 < Q < 3 GeV and (b) 4 < Q < 5 GeV. Sohd and 
dashed lines are for Re and Roir, respectively. 
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FIG. 15. The ratios in Eqs. (^) and ( |6l| ) as a function of Qt at y = for the LHC energy 
VS = U TeV and mass intervals (a) 2 < Q < 3 GeV and (b) 4 < Q < 5 GeV. Sohd and dashed 
hnes are for Re and Roir, respectively. 
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FIG. 16. Ratio Rg in Eq. (|6^ ) as a function of Qt at y = at the Tevatron energy (sohd), the 
LHC energy (dashed), and the RHIC proton-proton energy (dotted) for virtual photon invariant 
masses (a) Q = 2 GeV and (h) Q = 5 GeV. 
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